Experimental Parameters for a Reactor Antineutrino Experiment at Very 

Short Baselines 



K.M. Heeger, and M.N. Tobin 

Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

B.R. Littlejohif] 

Physics Department, University of Cincinnati, Cincinnati, OH 45221, USA 

H.P. Mumm 

National Institute of Standards and Technology, Gaithersburg, MD 20899, USA 
(Dated: December 11, 2012) 

Reactor antineutrinos are used to study neutrino oscillation, search for signatures of non-standard 
neutrino interactions, and to monitor reactor operation for safeguard applications. The flux and 
energy spectrum of reactor antineutrinos can be predicted from the decays of the nuclear fission 
products. A comparison of recent reactor calculations with past measurements at baselines of 
10-100 m suggests a 5.7% deficit. Precision measurements of reactor antineutrinos at very short 
baselines 0(1-10 m) can be used to probe this anomaly and search for possible oscillations into sterile 
neutrino species. This paper studies the experimental requirements for a new reactor antineutrino 
measurement at very short baselines and calculates the sensitivity of various scenarios. We conclude 
that an experiment at a typical research reactor provides 5(j discovery potential for the favored 
oscillation parameter space with 3 years of data collection. 



I. INTRODUCTION 



Recent calculations of the predicted Z/g flux from reactors compared to past measurements at baselines 
between 10-100 m have revealed an apparent deficit of about 5.7% [T]. This discrepancy, the "reactor 
anomaly" , can be interpreted as a sign of new physics or point to an issue with the reactor flux calculations 
[2]. Independent calculations have verified the reactor Ve flux predictions [3] but subsequent calculations 
based on a full ah initio prediction using newly available nuclear data |4j reduce the size of the discrepancy. 
It has been suggested that the ^^reactor anomaly'^ may be the signature of additional sterile neutrino states 
with mass splittings of the order of ^ IcV^ and oscillation lengths of 0(3 m) [S]. Other anomalies in neutrino 
physics including the observation of apparent Ve and appearance at similar mass-squared splittings in 
accelerator experiments [51 [7] , deficits in observed events from high- intensity i^e sources used to calibrate 
solar neutrino detectors 0, and the preference for more than three relativistic species in astrophysical 
surveys [S] add to the puzzle. 

Additional data will soon be provided by the currently-operating km-scale reactor experiments, Daya 
Bay, Double Chooz, and RENO [TU1 - [T^ . While highly precise, these experiments cannot resolve the short 
oscillation lengths associated with eV^ mass splittings; at these distances and with a finite detector energy 
resolution, the oscillation effect from potential sterile states averages to yield an effective rate deficit. 
Moreover, these measurements will eventually be limited by the understanding of the interference of multiple 
reactor cores, the presence of oscillation effects from other mass-squared splittings, and by the inability to 
take background data free of reactor Vg. A new experiment at very short (<10 m) baselines with a single 
core in a controlled research environment where backgrounds can be measured independently is needed to 
fully disentangle reactor flux and spectrum prediction uncertainties from sterile neutrino oscillation effects, 
or other signs of new physics. 
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The paper is organized as follows: Section [IT] summarizes both the experimental detector and reactor 
parameters under consideration. Section [lll| introduces the characteristic signature of neutrino oscillations 
and defines the analysis used to calculate experimental sensitivity. A nominal, generic reactor-detector 
arrangement referred to as the "default experiment" is defined for systematic studies and comparison. 
Vl 



Sections IV 



VI 



and VII present a discussion of reactor, facility, background, and detector parameters 
and examine the impact on the overall sensitivity of the experiment. Section [VIII| summarizes the sensitivity 
and discovery potential of a new reactor experiment and discusses its optimization. 



II. EXPERIMENTAL PARAMETERS 



A number of experiments have been proposed to address the reactor anomaly by making a precise mea- 
surement of the reactor flux and spectrum at very short baselines [T5Ul6j . The proposed experimental 
configurations are site-specific and refiect the infrastructure and logistical constraints at each site. This 
work studies generic experimental parameters that determine the ultimate sensitivity of a reactor experi- 
ment at very short baselines, and helps guide the optimization of a new experiment aimed at a definitive 
sterile neutrino oscillation search. 



A. Reactor Parameters 



Antineutrinos from the reactor core are used as the flavor-pure source. The relevant reactor reactor 
core parameters are the following: 

I. Reactor power: Each nuclear fission initiates the release of a known amount of energy along with 
several antineutrinos. The reactor flux is proportional to the thermal power output of the reactor 
modulo corrections for the isotopic fuel composition. High operating power combined with high reactor 
up-time maximize the event statistics. Most commercial reactors, but not all research reactors, operate 
close to their licensed power. 

II. Fuel type: Most nuclear reactors utilize uranium-based fuel containing a mixture of ^ssy g^jj^ ^35^^ 
with the latter providing the majority of total fissions and thermal power. The two most common classes of 
reactor fuels, commercial and highly-enriched, differ primarily in the total percentage enrichment of ^'^^U. 
Highly-enriched uranium (HEU) fuel, commonly utilized in research reactors, contains upwards of 90% ^■^^U, 
with fission fractions dominated by ^ssy g^jj points in the fuel cycle. Commercial fuel generally contains 
less than 6% ^'^^U. While ^35^ fission comprises the majority of all fissions in these cores, a significant fraction 
(30-50%, depending on the fuel burnup) is contributed by ^38^ and by ^'^^Pu and ^^^Pu accumulated during 
the fuel cycle |17j . These isotopes produce different proportions of fission products, resulting ultimately in 
a difference in the total number, energy spectrum, and spectral uncertainty of produced Ve- 

HI. Duty cycle: The duty cycle is given by the ratio of the power cycle to the subsequent shutdown 
periods for refueling or maintenance. The duty cycle influences the statistical power of an experiment 
in a straightforward way. Therefore, throughout this paper, we present results only in terms of detector 
livetime. Importantly, reactor down-time can be used for determining the shape and position distribution 
of backgrounds [TS] . 

IV. Core dimensions: Antineutrinos are produced throughout the active reactor core and emitted isotrop- 
ically. The core's finite dimensions cause a spread in the neutrino path lengths between the reactor core 
and detector that washes out the observable oscillation in the detector. 

As examples, we use parameters from three research reactors in the US, the High Flux Isotope Reactor 
(HFIR) at Oak Ridge National Laboratory (ORNL) the National Bureau of Standards Reactor (NBSR) 
at National Institute of Standards and Technology (NIST) |20j, and the Advanced Test Reactor (ATR) 
at Idaho National Laboratory (INL) Jii. The Institut-Laue Langevin (ILL) reactor in France and the 
commercial San Onofre Nuclear Generating Station (SONGS) in California are included for reference [^123) . 
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B. Facility Parameters Parameters 

The reactor facility provides constraints on the location of the detector, its size, and distance to the core. 
The following parameters determine the experiment's design and sensitivity: 

I. Experimental space and detector volume: Detector location, geometry, and dimensions are limited 
by the available space near the reactor core. Typically, biological shielding limits the distance of closest 
approach to the reactor while the cross-sectional area of the detector and the maximum radial distance are 
constrained by the layout of the reactor facility, floor-loading limits, and the availability of floor space not 
occupied by other experiments, reactor equipment, or detector shielding. To maximize event statistics, a 
large cross-sectional area at close radial distances is desirable. 

II. Detector distance to core: The closest distance to the reactor determines the overall I'e flux seen by the 
detector as well as the oscillation wavelength probed. For the following discussion we define this distance of 
closest approach, r„ji„, to be from the center of the reactor core to the closest point in the active detector 
target. 

Throughout this discussion we explore a range of detector sizes and closest distances reflecting the ex- 
pected space available at US facilities. 

C. Detector Parameters 

Intrinsic detector characteristics determine the observed event distributions: 

I. Fiducial volume and target mass: The observed number of Jy^ interactions in the detector scales with 
the number of target protons in the detector and sets the statistical precision of the experiment. 

II. Detection efficiency: The detection efficiency directly scales the detected number of events. For fixed 
signal-to-background, S/B, changes in the detection efficiency have the same effect as increasing the detector 
volume or target mass. 

III. Position and energy resolution: An unambiguous demonstration of neutrino oscillation requires the 
observation of oscillations in energy and distance, the characteristic L/ E dependence. Good position and 
energy resolution are necessary to maximize the experimental sensitivity. 

D. Backgrounds 

A reactor experiment at very short baseline requires the operation of a detector on the surface under min- 
imal overburden. Muon and cosmic ray-induced backgrounds as well as backgrounds from the reactor are 
important considerations. Both the magnitude and shape of backgrounds are expected to be site-specific. 
Reactor-correlated backgrounds are affected by the reactor-detector distance, the reactor power, and the 
amount of shielding material between the reactor and detector. Cosmic background rates are dependent on 
the total overburden provided by the facility and the mass distribution near the detector location. Back- 
grounds can be mitigated with active or passive shielding, as well as with background rejection techniques, 
such as pulse-shape discrimination or detector segmentation. Small but irreducible backgrounds similar in 
energy spectrum to the signal may also be provided by spent nuclear fuel nearby the reactor [24| . Detailed 
background studies are not presented for any particular reactor sites in this paper; instead, reasonable 
background models are constructed based on experience in previous experiments [221 [^5H?7] . Specifically, 
we consider the following parameters: 

I. Signal-to-background ratio (S/B): The overall S/B ratio is a measure of the total magnitude of time- 
coincident backgrounds and accidental background in the delayed coincidence window of the inverse beta- 
decay reaction at all energies and positions. 

II. Background spectral shape: The backgrounds' energy dependence will determine the measured spectral 
shape of the observed events. Energy-dependent background subtraction and/or fitting will be important 
for the analysis of the energy-dependent oscillation signature. 

III. Background position distribution: In the vicinity of a reactor significant spatial variations of back- 
grounds are expected. In particular, fast neutrons can scatter from surrounding materials, building struc- 
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tures, and even other experiments. Measurement of the observed event rate as a function of position through 
detector movement or segmentation wiU be critical for understanding local background variations on the 
meter-scale. 



III. EXPERIMENTAL SIGNATURE 



A. Oscillations of Reactor Antineutrinos 



Antineutrinos from reactors are produced as flavor-pure V,, in the decays of the neutron-rich fission 
products in the reactor fuel. More than 99.9% of all emitted from commercial reactors are produced 
within the decay chains of four isotopes, ^^^U, ^^^U, ^^^Pu, and ^"^^Pu. The thermal heat released in the 
nuclear decays is proportional to the number of emitted and is thus a measure of the flux of expected 
antineutrinos. The spectrum of reactor antineutrinos detected via inverse beta decay has a mean energy of 
about 4 MeV and extends up to roughly 10 MeV. 

Neutrinos and antineutrinos are produced as a linear combination of mass eigenstates and their fiavor is 
associated with the accompanying lepton. Due the difference in the mass eigenstates the flavor of observed 
neutrinos oscillates as a function of baseline and energy. For the three active neutrino states the neutrino 
mixing parameters are well measured in atmospheric, solar, reactor, and accelerator based experiments. Re- 
actor Ve. disappearance over baselines of 1-2 km and 180 km has been observed. The oscillation probability 
can be parameterized in terms of the mass splitting Am|^ and the mixing angle 0^ between the ith and jth 
mass eigenstate. Additional sterile neutrino mass states with Am^ ~ 1 eV^ beyond the 3 active neutrinos 
would yield an oscillation effect over meter-long baselines with survival probability described by 



P,^r{E, L)^l- sin^ 26,, sin^ [ l-27AmL|£ r\ \ 



with oscillation parameters Am|]^ and sin 29, 




FIG. 1: Unoscillated (left) versus oscillated (right) detected spectra at (Am^ = 1.8 eV^, sm^29^^ = 0.5) for 
a reactor Ue detector with realistic parameters as described in Table |T] Exaggerated oscillation values are chosen 
for illustrative purposes. The change in spectral shape from baseline to baseline is a key signature of neutrino 
oscillations. 

Figure [I] illustrates the oscillation effect in baseline and energy for (Am2=1.8 eV^ sin^ 26*^6=0. 5). 
The characteristic L /E oscillation is pictured in Figure for the sterile neutrino oscillation parameters 
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(1.8 eV^jO.l) preferred in global fits [5] and the default experimental arrangement described below. 
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FIG. 2: Ratio of the oscillated to unoscillated spectrum as a function of L/E for (Am^ = 1.8 eV^ sin 2^ee=0.1) 
and the nominal experimental arrangement described in |III B[ Error bars are purely statistical. 



B. A Nominal Reactor Antineutrino Experiment at Very Short Baselines 

Reactor neutrino experiments typically utilize the inverse beta decay reaction V,, +p — )• e+ +rt (IBD) with 
a threshold of roughly 1.8 MeV to measure the flux and energy spectrum of reactor V^- Liquid scintillator 
detectors provide a proton-rich target with high detection efficiency and good energy resolution. Rejection 
of backgrounds is achieved by means of the delayed coincidence signature and correlating cosmic rays with 
muons-induced events. Energy resolutions of around 8% and total detection efficiencies of 75-85% have been 
obtained in recent large-scale underground experiments (TUl E] . For smaller detectors developed for early 
short-baseline experiments and more recent reactor monitoring purposes, energy resolutions of 10-20% 
and efficiencies of 10-50% have been reported p51[25ti5U] . For the studies presented in this paper a nominal 
experimental arrangement of C'(m'^)-sized detector at distances of 0(10 m) from an average research reactor 
is considered. The specific parameters assumed for the nominal experiment are listed in Table |lj 



C. Experimental Sensitivity and Discovery Potential 

The sensitivity of a reactor experiment to neutrino oscillations is evaluated by comparing the detected 
energy spectrum to the expected one in the absence of neutrino oscillations. A radially extended detector 
with position resolution allows a comparison as a function of baseline. A test is used to test the hypothesis 
of no-oscillation and for parameter estimation in Am|j^ and sin^20ee- Observed V^, events are binned in 
energy E with index i and in distance between the core and the point of detection, \L\, along index j. The 
expected unoscillated number of events per bin, Tij, is given by 

mr).^/^fl^.., (2) 

with Np as the number of target protons, the detection efficiency e, energy E, point vector f* between the 
core center and production point in core, vector L between the core center and detection point in the 
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Parameter 


Value 


Comment 


Reference 




Power 


on Aff\A?" 
zU IvlW 


INib i-iike 


|oi 1 


Reactor 


Shape 


cylindrical 


iNib i-iike 


rm 




oize 


0.5 m radius, half- height 


IN io i -ilKe 


Iqi 1 




ruei 


PTTTTT 

ritii\j 


Research reactor fuel type 


IQI NQOI 




Dimensions 


i X i X o m 


3 meters of available baseline 




Detector 


Efficiency 


OU /O 


111 1 dll^c (Jl OlDlj cXpb. llU-tJU/0) 


rj^ F7n\ 

IZtJI IZiOl |ZiC7l 


Proton density 


6.39x10^**^ 


From Daya Bay GdLS 


m 




Position resolution 


15 cm 


Daya Bay-like 


mi 




Energy resolution 


10%/Ve 


Daya Bay-like 


m 




Run Time 


1 year live-time 






Other 


Closest distance 


4 m 


NIST-like 






S:B ratio 


1:1 


In range of SBL exps. (1-25) 






Background shape 


l/E^ 


Similar to SBL experiments 


[smsiiiz] 



TABLE I: Nominal experimental parameters used for the sensitivity calculations presented in this paper. 



detector, energy spectrum S{E), and inverse beta cross-section cr{E). Mij is the expected number Tij 
with backgrounds added in the presence sterile neutrino oscillations as described by Equation IT] Both Mij 
and Tij are subject to gaussian position and energy resolution smearing according to the values given in 
Tabled) The function is defined as: 



[M,j -{a + al+ apT.j - (1 + ab)Bij] ^ 
Tij + {ab2bBijY 




The sums over 17 prompt energy and 19 position bins in the range of [0.8,7.6] MeV and [3.2, 7.4] 
m, with bin withs of 0.4 MeV and 0.2 m respectively. The bin widths are comparable to their respective 
modeled detector resolutions (10%/-\/i?, 0.15 m). The sum is minimized with respect to 9ee^ Am|]^ and 
to the nuisance parameters {a, a^.^ as described in 36J. The parameter a allows the signal 

normalization to vary within the bounds of its associated uncertainty a to account for uncertainties in the 
absolute reactor V,, normalization and absolute detection efficiency. The 100% error in a ensures a floating 
overall normalization, meaning sensitivity to oscillation is only given by spectral distortions. 

The parameters a* account for the the uncertainty in the reactor v^. spectrum from reactor flux predictions 
and from previous experimental measurements, as well as detector systematics uncorrelated between energy 
bins. These parameters allow position bins at one energy to fluctuate together, independently of position 
bins at any other energy. These correlated fluctuations are limited by cr*, which vary as a function of energy 
as described in j3]. Modelling of these uncorrelated errors in energy is of particular importance, as they 
limit the power of a pure energy-based oscillation analysis. This will be discussed further in Section [VlIIB[ 

The position spectrum parameters al allow correlated fluctuations with position, rather than energy, in 
order to incorporate the effects of relative efficiency differences and uncertainties between position bins. In 
contrast to the energy spectrum uncertainties, relative efficiency uncertainties between position bins, given 
a value CTr=0.5% for all bins, should be smaller, as they are easier to characterize via detector simulation 
and calibration. 

Backgrounds Bij are estimated with a flat position dependence and an energy spectrum that falls with 
energy as generally mirroring the background shape reported by some previous short-baseline exper- 

iments [m . The background normalization is allowed to fluctuate similarly to the signal normaliza- 

tion by incorporating the background nuisance parameter a?,, with an associated systematic uncertainty of 

(Tfc = 10%. 

Uncorrelated uncertainties in the background energy and position spectrum should be very well- 
characterized for a successful experiment. The sources of uncertainties in position spectrum and in cn- 
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ergy spectrum for backgrounds are not likely to be mostly decoupled, as is the case for the signal position 
and energy spectrum uncertainties. Without the specifics of these correlations, the most conservative way 
to incorporate such an uncertainty into the analysis is to provide an additional independent nuisance 
parameter for every bin with an associated uncertainty ab2b that reflects the precision of any background 
spectral measurements. Conservative consideration of these uncertainties can be more simply achieved by 
adding the effect of ab2b to the denominator of the x^- For this study, (T(,2& is given a default value of 0.5%. 

The 3a and 5a discovery potential for neutrino oscillations is calculated following the prescription in |36) . 
The resulting discovery contours are shown in Figure [3] The nominal reactor experiment considered in 
Table |l] is capable of excluding a large fraction of the currently preferred parameter space to better than 3a 
with one year of live-time and at 5a C.L. with 3 years of data taking. Sections |IV[ \V\ and |VII| investigate 
the impact of the reactor, facility, and detector parameters on the sensitivity of the experiment. The 
experimental parameters listed in Table |l] will be used throughout this paper as reference for the nominal 
experimental arrangement. 
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FIG. 3: Discovery potential to neutrino oscillations at 3- and 5a for the nominal experimental arrangement described 
by Table [I] The 3cr contour corresponds to one year of detector live-time while the 5a contour is shown for three 
years of detector live-time. Also pictured are the best-fit parameter spaces for the reactor antineutrino anomaly [5] 
and for the 3-1-1 global fit to all relevant accelerator, source, and reactor data given in |37| . 



IV. REACTOR PARAMETERS 



A. Reactor Power and Operations 



The flux of z7e emitted from reactors is directly related to its thermal power. The total nominal power 
capacities for several research and commercial reactor sites are shown in Figure |4] At 3GWt;i, typical 
commercial power stations such as the San Onofre Nuclear Power Generating Station (SONGS) are 
about an order of magnitude more powerful than research reactors. Research reactors in the US have a wide 
range of thermal powers up to a maximum of 250 MW at the Advanced Test Reactor [33] • The variation of 
the experimental sensitivity of the baseline experimental conflguration with the power of several research 
reactors (or Ve flux) is shown in Figure |4j The increase in event statistics with thermal power uniformly 
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increases an experiments' sensitivity for all values of Am^ 
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FIG. 4: Left: Comparison of nominal and operational thermal power output for selected reactor facilities. The 
nominal licensed power output (dashed) is roughly equal to the operational power (solid) for all facilities except the 
Advanced Test Reactor. Right: Sensitivity of the nominal experimental configuration for various power levels and 1 
year of detector livetime. 

Nuclear reactors are periodically shut down for refueling and maintenance. As a result, the time-averaged 
thermal output of a nuclear reactor will be less than its maximum thermal output over that time period 
because of periods of below-capacity operation and reactor-off periods. For example, at ATR the standard 
operating thermal output of 110 MW is significantly below its maximum licensed thermal output of 250 MW. 
See Figure [4] All nuclear facilities undergo reactor-off periods in which reactor refueling and maintenance 
takes place. Table [iTl summarizes the length and frequency of reactor-off periods at several reactor facilities 
as well as their reactor-on and -off time. 



Reactor 


Power 
(MWih) 


Baselines 
(m) 


Reactor On 
(Days) 


Reactor Off 
(Days) 


Down- Time 


Ref. 


NIST 


20 


4-13 


42 


10 


-32% 


m 


HFIR 


85 


6-8 


24 


18 


-50% 


m 


ATR 


250 (hcensed) 
110 (operational) 


7-8 (restricted) 
12-20 (full access) 


48-56 


14-21 


-27% 




ILL 


58 


7-9 


50 


41 


-45% 


m 


SONGS 


3438 


24 


639 


60 


8.6% 


[MIISS] 



TABLE II: Summary of reactor powers, accessible baselines, fuel cycles, and reactor-off times at various reactor 
facilities. At ATR, detector access is limited for baselines <12m. SONGS data is from past operation; this facility 
is currently shut down. The down-time includes estimates for seasonal shutdowns and maintenance periods. 

Reactor-off time and operation at reduced power directly reduce the total annual z7g event statistics 
yielding the same net effect as a lower thermal power capacity. The experiment's sensitivity decreases with 
lower power as illustrated in Figure [4j However, reactor-off time provides an opportunity to measure the 
rate and energy distribution of backgrounds. A detailed understanding of the spatial and energy distribution 
of backgrounds is critical for a precision experiment at short distances from the reactor. This topic will be 
discussed in greater detail in Section |VI[ The optimum ratio of background measurement time relative to 

signal time will depend on the total signal statistics, the dominant uncertainties, as well as the signal- 
to-background ratio. This optimization should be carried out as part of a detailed design process. For the 
research reactors in Table |TT1 the I7g signal time will be between 50-70% of calendar time. 



9 



B. Reactor Fuel 



Commercial nuclear power stations use conventional nuclear fuel comprised of a mixture of U and Pu 
isotopes while some research facilities operate with highly-enriched uranium (HEU). The four isotopes ^■^^U, 
238u^ 239p^^ g^j^jj 2'*iPu produ ce > 99.9% of all I'eproduced in a reactor. In HEU cores nearly all fissions are 
accounted for by ^'^^U. Table III gives the main fission isotopes in the two reactor fuels and their relative 
contributions to the total fission rate. 



Fuel Isotope 


Time- Averaged Fission Fraction 


Conventional Fuel 


HEU fuel 


235 


0.59 


>0.99 


238 y 


0.07 


<0.01 


239p^ 


0.29 


<0.01 


2«PU 


0.05 


<0.01 



TABLE III: Approximate time-averaged fuel compositions for various reactor cores. Fractions for conventional |17| 
and HEU reactors [IT], respectively. 



The Ve spectrum and rate per fission is different for each of these isotopes |42) and the flux and spectrum of 
reactor V,, are the sum for all isotopes in the reactor core. The time-averaged detected V,, spectra from HEU 
and conventional fuels are compared in Figure [sj The integrated flux differs by roughly 8%, and the time- 
averaged spectral differences are 10% or less. The time evolution of the isotopic fuel composition creates 
a time-dependent spectral shape. For the spectral range shown in Figure [5] the typical fission fractions 
for the average fuel were evaluated and compared to the upper and lower fractions at the beginning and 
end of fuel cycle [T7]. Using the spectral shapes from [32] together with the isotope fraction the combined 
spectrum can be calculated. When convoluted with the standard IBD cross section we obtain the detected Vg 
spectrum. The differences in the energy spectrum only have a small impact on the experiment's sensitivity 
to short-baseline neutrino oscillation as shown in Figure |6] 




FIG. 5: Left: Detected reactor spectra from conventional commercial reactor and HEU reactor fuel (top panel), 
and their spectral differences (bottom panel). The band indicates the change in the spectral shape over one fuel 
cycle. Right: Fractional change in spectrum over one fuel cycle for conventional and HEU reactors. 

In addition to differences in the integrated flux and spectral shape, the uncertainties spectral shape vary 
between the dominate fission isotopes. These arise from a combination of the statistical and measurement 
uncertainties of fission isotope beta spectra |22j and uncertainties in the conversion of electron spectra 
to corresponding Vg. spectra [311321 133] ■ Spectral uncertainties for ^ss^j j-^nge from 1.8-3.2% in the range 
from 2-6 MeV, while for ^sspu they increase to 1.9%-5.7% and 2.5%-5.0% for 24ipu [35]. Because of these 
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FIG. 6: Sensitivities of the nominal reactor experiment for conventional and HEU fuel reactor cores. The use of 
HEU fuel leads to minor improvements in the overall sensitivity. 



considerations, spectral uncertainties will be lower for HEU fuel than for conventional fuel. The change in 
sensitivity resulting from this difference in spectral uncertainty for the two fuel types is shown in Figure |6] 
HEU fuel provides a minor improvement in sensitivity as a result of the lowered spectral uncertainties. The 
benefits of lower spectral uncertainties will be amplified in the case of detectors with limited baseline ranges 
or position resolution, as will be discussed in Section [VHI B[ 



C. Reactor Core Size and Dimensions 



Reactor cores come in a wide variety of shapes and sizes. The production in a reactor follows closely 
the distribution of the fuel assemblies. For the purpose of a neutrino oscillation experiment at very short 
baselines the main difference between commercial power reactors and research reactors is in their size. 
Conventional pressurized water power reactors are typically 3-4 m in height and diameter, while research 
reactors tend to be smaller in size with dimensions as compact as ^0.5 m or less. Figure [7] provides an 
illustration of the variety of shapes and sizes of various research and conventional reactor cores. 

For neutrino oscillation searches it is required that the size of the reactor core and thus the spread in 
neutrino path lengths to be less than the oscillation wavelength to avoid a wash-out of the oscillation signal. 
Figure [7| shows the spread in path lengths between the finite-sized reactor cores and a point-like detector 
horizontally displaced from the vertical midpoint of each reactor at a distance of r = 10 m. The distance, 
r, is defined between the geometric center of the core and the point-like detector. A distance of 0(10 m) 
represents a typical distance for very short baseline reactor experiment. We assume that neutrinos are 
produced and emitted uniformly throughout the core region. We define the path length, I, as the distance 
between the points of production and observation. Convolving the path length distribution with l/R^, 
determines the expected spatial distribution for the relative probability of Vg interactions in the detector. 
We note that a detector with finite position resolution or segmentation adds additional smearing to the 
observed path length distribution and the observed oscillation effect. 

For finite-sized reactor cores of any geometry we can calculate the average path length, I and the RMS 
of the path lengths as 



(4) 
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FIG. 7: Left: Approximate size and shape of fuel distribution in five reactors including the US research reactors 
ATR, NIST, and HFIR. SONGS is a conventional power reactor, for comparison. Right: Relative event probability 
versus neutrino path lengths for several finite-sized reactor cores and a point-like detector horizontally displaced 
from the vertical midpoint of each reactor at a distance of r = 10 m. 
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Figure [S] shows the average neutrino path length spread for different core geometries as a function of 
distance from the reactor. At d > 5m the RMS spread in path lengths approaches ^0.5-0.6 m for a reactor 
of 1 m height and diameter. The variation due to the core shape is significantly smaller than the total 
magnitude of spread. For all practical considerations and for realistic distances of d > 4 m from the reactor 
the shape of the core and fuel distribution is only of secondary consideration. For highly asymmetric cores 
such as long, cylindrical arrangements choosing the orientation of the detector with respect to the symmetry 
axis of the core can be used to reduce the spread in neutrino path lengths. Figure |9] shows the discovery 
potential for the nominal reactor experiment for varying dimensions of the reactor core. The sensitivity to 
higher Am^ values is lost as the core width is increased. The overall core dimensions dominates the spread 
in neutrino path lengths while the reactor shape only plays a secondary role. 
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FIG. 8: Left: Geometries, cross-sections, and path lengths for a cylindrical reactor core and point-like detec- 
tor. Right: Corresponding path length spreads for a 1-m cylindrical core. The differences in path length spreads 
contributed by the circular horizontal dimensions and the rectangular vertical dimensions of the core are a small 
correction on top of the overall magnitude of the spread, which is defined by the overall core size. 
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FIG. 9; Variations in the discovery potential of the nominal reactor experiment with a half-height core of various 
radii (left) and core geometries (right). The overall spread in neutrino path lengths is dominated by the overall 
dimensions of the core. 



V. FACILITY PARAMETERS 



A. Experimental Area and Detector Volume 

A key signature of neutrino oscillation is the distortion of the measured spectrum with energy and the 
variation of the energy spectrum with baseline. The baseline dependence of the energy spectrum can 
only be observed in a detector whose active length comprises more than a small portion of one oscillation 
period. Ideally, one would map out the oscillation over an entire wavelength or more. The availability of 
experimental space in a reactor facility limits the accessible baselines as well as the total active detector 
volume. The accessible baselines listed in Tabic [Til provide a distance range of up to 9 m at one particular 
facility. This allows the placement of point-like detectors in various locations, or the construction of an 
extended detector in the radial direction. For any significantly extended detector the change in event rate 
due to the 1/R^ law has to be taken into account. 

An increase in detector length has two effects: It increases the total overall event statistics of the experi- 
ment and it samples a larger fraction of the possible oscillation wavelength. Both the maximum achievable 
sensitivity and the overall range of accessible Am^ values improve with an increase in detector length. In 
particular, sensitivity to lower values of Ato^ is improved with a longer detector as oscillations with longer 
wavelength are sampled. Figure [TO] shows the change in sensitivity of an experiment with the default char- 



acteristics for three different detectors of 1, 3, and 5 m lengths. Figure 10 shows on one side the overall 
sensitivity of the three detectors with the different event statistics due the various detector lengths and on 
the other hand the same experimental arrangement normalized to the statistics of the default experiment. 
The normalized case illustrates the impact of the detector length or baselines while the un-normalized 
situation highlights the decrease in event rate with l/B?. 

The detector cross-section is typically limited by space constraints inside the containment building near 
the reactor core. As the overall target mass and thus statistics at each baseline scale with the cross-sectional 
area of the detector, one can consider scaling the cross-sectional area of the detector as a function of baseline 
to counteract the effect of the 1/i?^ reduction in flux. 
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FIG. 10: Variation in sensitivity of the default experiment with detector length. Left: Event statistics is normalized 
to the default experiment in each case. Both the maximum achievable sensitivity and the range of accessible Am^ are 
increased with increased detector length. Right: Un-normalized case demonstrates the additional effect of increased 
statistics with larger detector volume. 



B. Reactor-Detector Distance 



For extended core and detector geometries, the distance between the reactor and the detector is not 
uniquely defined. We therefore choose the reactor-to-detector distance, r, as the distance between the 
center of the reactor core and the closest point to the reactor in the active detector region while the detector 
length, d, describes the radial length of the active detector volume. The reactor-to-detector distance, r, at a 
minimum, is comprised of the extent of the reactor itself and the thickness of the containment and shielding 
structures surrounding the core. In practice some passive shielding will be required to operate a V,, detector 
in the vicinity of a reactor. The accessible baselines for various reactor facilities listed in Table |ll] take 
0.5 m additional space into consideration for this purpose. The closest accessible baseline is site specific 
and referred to as r-mm • As is evident in Table [Ilj research reactors provide an opportunity for improved 
measurements of the reactor X/g flux and spectrum at the shortest baselines to date. Facilities such as NIST 
may provide access to baselines as short as 4 m or less. 

In maximizing event statistics and experimental sensitivity, the design of an experiment is a trade-off 
between distance, spread in neutrino path length, and relative reactor power. Figure [TT] shows the spread in 
neutrino path lengths from finite-sized reactor cores as seen by a point-like detector at the closest accessible 
baseline rmm- The spread in neutrino path lengths reflects the size of the reactor cores while the area of 
each distribution indicates the relative reactor power scaled by . 

The right panel of Figure [TT] shows the sensitivity of the default experiment for Tmin = 4, 6, and 12 m 
respectively. To emphasize the effect of varying rmin, the normalization is adjusted to maintain the same 
total statistical sensitivity in each case. The closest accessible baseline, r„iim impacts the high-Am^ sensi- 
tivity of the experiment. This is due to the fact that finite energy resolution tends to wash out the observed 
oscillation with distance from the reactor core. 

The effect of rmin on the experiment's sensitivity to various ranges of Am^ is dependent on the total 
detector length d. Longer detectors help improve the sensitivity to lower values of Am^ but do not compen- 
sate completely in case of larger r^m- This is illustrated in Figure [T2] for a 3 m- long detector placed at 4 m 
and 12 m respectively from the reactor core. Figure [T2| illustrates how a detector of chosen length samples 
different fractions of the oscillation periods and oscillation amplitudes for various Am^. For larger Am^, 
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FIG. 11: Left: Relative probability of interactions in a point-like detector at the closest possible distance from 
the respective reactor cores. The spread in neutrino path lengths reflects the size and shape of the reactor cores. 
The peak areas give the relative reactor power scaled by 1/B?. Right: Sensitivity of the default configuration for 
closest accessible baselines of 4, 6, and 12 m respectively. To emphasize only the effect of baseline, the normalization 
is adjusted to provide equal total event statistics. 



the best achievable sensitivity at large rmin is lower even with longer detectors because the amplitude of 
the oscillation effect is diminished. For low Am^, only a small portion of the oscillation period fits inside 
the detector length at both distances, indicating the need to sample a much longer range of baselines. For 
the favored mass splitting the oscillation period is sampled well at both chosen baselines. 

If the oscillation length were known a priori, a detector of suitable length could be placed between 
the first oscillation minimum and subsequent maximum to maximize the observed oscillation difference or 
around the oscillation minimum to observe the turning point. Ideally, such a detector would be movable to 
measure both the difference between the oscillation maximum and minimum and the turning point around 
one of the oscillation extreme. A movable, extended detector can also help mitigate possible backgrounds 
or systematic effects that can mimic this oscillatory signature. In the case of longer wavelengths and given 
the facility constraints at research reactors, multiple, radially extended detectors may be necessary to allow 
for a comprehensive search and discovery of neutrino oscillations with unknown /S.m? . 
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FIG. 12: Oscillation of 4MeV for Am^ of 6.0, 1.8, and 0.15 eV^ as a function of baseline. The shaded regions 
indicate the position and radial dimension of the detectors for the default experimental configuration. The dimensions 
of the reactor core are shown in solid color. Oscillation curves include position smearing resulting from position 
resolution and from the finite reactor core size. Boxed areas indicate the range of oscillated/unoscillated oscillation 
amphtudes covered by each case. For the large Am^ value (5 eV^), the 3m detector length encapsulates more 
than a full period, and the damping of the oscillation amplitude due to finite detector resolution reduces the visible 
oscillation at farther distances. For medium Am^ values (1.8 eV^), the 3m detector length allows for significant 
observation of oscillation a both distances. For very small Am'^ (0.15 eV^), the variation in oscillation from detector 
font to back is larger at 12 m than at 4 m, although both differences are small, as the 3 m detector length encapsulates 
only a small portion of the total period. Two detectors at different baselines may be needed to probe regions of 
small Am^. 



VI. BACKGROUNDS 



The overall magnitude and spectral shape of backgrounds in a short-baseline reactor neutrino experiment 
are determined in large part by the detector's surroundings, including the distance to the core, nearby spent 
fuel repositories, neutron backgrounds from the reactor and nearby experiments, and cosmic ray induced 
backgrounds. Also important are the cleanliness of the detector components, as well as the use of any 
background-reducing detection detection techniques. The exact shape and magnitude of the backgrounds 
may vary widely and will only be known by conducting site-specific background surveys, detector material 
surveys, and detector simulations. 

In place of detailed background surveys, the sensitivity of the default experiment to varying background 
conditions was investigated by varying the magnitude and shape of the input background spectrum in the 
calculation. Signal-to-background ratios of 0.1, 1.0, and 10 were considered. Changes in sensitivity resulting 
from variation of S:B ratio of the default experiment can be seen in Figure [TSj The overall background rate 
clearly has a significant impact on experimental sensitivity at all values of Am^. 

As defined in the given in Equation [s} backgrounds are assigned a nuisance parameter that accounts 
for the uncertainty in the overall background normalization. Wide variation of this parameter's associated 
systematic uncertainty has negligible effect on the sensitivity of the default experiment. 

Background spectral shapes were also varied in ad diti on to overall background normalization. Three dif- 



ferent background spectral shapes, pictured in Figure 14 were used: the default XjE"^ shape seen in previous 



short-baseline experiments l22ll25U27j . a flat distribution commonly associated in deep-underground reactor 
Vf, experiments with fast neutron backgrounds, and a distribution identical to the signal distribution. The 



results of these variations are also shown in Figure 13 The overall background normalization appears to 



have a much larger effect on sensitivity than the exact spectral shape, within the assumptions of this study. 
The same is also found to be true for the shape of the background position distribution. 

As discussed in Section |III C[ uncertainties in the background energy and position spectrum shapes, 
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FIG. 13: Variation in sensitivity of the default experiment with overall signal-to-background (left) and shape (right). 
The default parameters used throughout the paper are a S:B ratio of 1:1 and a spectral shape. Backgrounds will 
play a critical role in the experimental sensitivity at a broad range of Am^ values, with the exact effect depending 
on the precise spectral shape as a function of energy and position. 




FIG. 14: Left: Example signal and background energy spectra. A SignahBackground ratio of 1:1 as used in the 
default experiment is shown. With a 1:1 ratio, the signal-shaped background spectrum is identical to the signal itself. 
A 5% uncertainty band indicates the effective range of uncorrelated variation allowed in the background spectrum 
by (Th2b- Right: Sensitivity of the default experiment to sterile neutrino oscillaions for various values of o"f,26. The 
significant variation in sensiivity with Oh2b underscores the importance of a precise knowledge of the background 
energy and position spectrum. 



without site-specific background surveys, conservatively accounted for with the addition of ah2b to the 
denominator of the . To determine the effect of this uncorrelated spectral uncertainty on overall sensitivity, 
0-626 was varied for the default experiment from 0.5% to 10.0%. The effect of this parameter is illustrated 
in Figure [M] each bin in energy and position is effectively allowed a free, uncorrelated fluctuation within a 
band determined by (Jb2b- For a real experiment, the magnitude of (7626 will be determined by measurements 
of background spectra during reactor-off periods, as well as measurements of low- and high-energy singles 
during reactor-on periods. Shifting the detector radially may also help map out and constrain the spatial 
variation of the backgrounds. The resultant change in sensitivity of the default experiment with variations 
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in cr{,2f) are also shown in Figure 14 The experiment's sensitivity is highly dependent on ab2b ioi all values of 
Am^. Precise knowledge of the background spectral shapes is clearly necessary for a definitive short-baseline 
oscillation experiment. 



VII. DETECTOR PARAMETERS 



A. Target Mass and Efficiency 



The total neutrino interaction rate scales proportionally with target mass and detection efficiency. Vari- 
ation of target mass and efficiency have an identical effect on the sensitivity as variations in reactor power. 
The target mass can be increased by increasing the cross-section of the detector, which is largely constrained 
by the available facility space, or by increasing the proton density of the target. Efficiency is increased by 
increasing the fraction of prompt and delayed signal energy deposited in the active scintillating region of 
the detector or by relaxing signal selection cuts. 



B. Detector Resolution and Analysis Binning 



Independent of a specific reactor-detector orientation and facility space constraints, the detector param- 
eters that define the sensitivity of an experiment are the detector's energy and position resolution. Good 
energy resolution is needed to both resolve the spectral distortions at a particular baseline and avoid washing 
out the oscillation at longer baselines, while position resolution allows for the observation of an oscillation 
as a function of distance. Previous short-baseline experiments have demonstrated the ability to attain sub- 
10% energy resolution [2H1I3D], while larger ^13 LS detectors have achieved resolutions of around 7-8% [55] . 
Acceptable position resolution can be obtained either through optical segmentation of the detector volume, 
or by using PMT charge topology information to reconstruct event positions in a larger one-zone detector. 
The former technique has been demonstrated in a number of previous short-baseline experiments, such as 
Palo Verde and Bugey [5S1I31], while the larger Gadolinium-loaded liquid scintillation detectors at Daya 
Bay, RENO, and Double Chooz have demonstrated successful position reconstruction with a resolution of 
30 cm or better pn[?7] . 

Figure [15] demonstrates the effect of variations in detector resolution for the default experiment. The 
fractional change in amplitude of the measured L/E oscillation pictured in Figure [2] is taken as the figure of 
merit. The oscillation signature is increasingly washed out as the position and energy resolution decrease. 
Attaining an energy resolution of better than 10% will not provide significant gains in sensitivity for this 
range of Am^. Position resolution approaching roughly 1/4 of the characteristic oscillation length at a 
particular Ato^ is required to maximize sensitivity. 

The effect of finite position and energy resolution can be approximated by varying the bin width in the 
analysis: the bin width provides an effective hard limit on the knowledge of exact event baselines and 
energies. The results of this approach are pictured in Figure [16] The average detections in each bin are 
given in Table IV to give a sense of the statistics provided by each binning. It is clear that sensitivity to 
high Am^ values is degraded as bin width is increased. 

Figure [16] also illustrates qualitatively the individual contributions of energy and position information to 
the total sensitivity. As the binning in energy or position, and thus the information available from energy 
or position, is reduced, sensitivity is reduced evenly, indicating that the two variables contribute roughly 
equally to the overall sensitivity of the default experiment. This is a feature of our choice of a default detector 
and the resolutions considered. Figure [16| also presents sensitivity for various binning choices of the same 
arrangement at 12 m closest distance. For this arrangement, sensitivity is reduced much more quickly when 
energy binning is reduced, indicating that energy information provides most of the experimental sensitivity, 
and that energy resolution is of paramount importance for longer-baseline detectors. 
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FIG. 15: Fractional change in maximum oscillation sensitivity with changes in energy (left) and position (right) 
resolution for the default experiment at either 4 m or 12 m closest distance with for Am^ values of 0.15, 1.8, or 5.0 
eV^ Fixed values of 10% for energy and 15 cm for position are used when varying position resolution and energy 
resolution, respectively. 
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FIG. 16: Variation in sensitivity to short-baseline oscillations with energy and position bin size for the default 
experiment at closest distances of 4 m (left) and 12 m (right). Statistics are normalized to the level of the default 
experiment, 250,000 Ve/year. Sensitivity is degraded, particularly at high Am^, as bin width is increased. One 
can also see the differing contributions to the total sensitivity from position and energy information, and how this 
contribution changes with differing detector closest distance. 



VIII. SUMMARY AND DISCUSSION 



A. Key Experimental Parameters for an Experiment at Very Short Baselines 

In this paper we have studied an optimized very-short-baseline oscillation experiment utilizing an HEU 
core with the highest possible power and smallest core size. To probe eV-scale mass-squared splittings, we 
consider a detector of at least 3 m at the shortest possible core distance with the highest possible efficiency, 
energy resolution of better than 10%, and position resolution of better than 0.2 m. Furthermore, a S/B 
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Number of bins 


Events 


Comment 


Energy 


Position 


per bin 


17 


21 


700 


Energy+Position Analysis, (0.4 MeV, 0.2 m) binning 


17 


6 


2450 


Energy+Position(Poor) Analysis, (0.4 MeV, 0.8 m) binning 


5 


21 


2380 


Energy(Poor)+Position Analysis, (1.6 MeV, 0.2 m) binning 


1 


21 


11900 


Position-only Analysis, 0.2 m binning 


17 


1 


14700 


Energy-only Analysis, 0.4 MeV binning 



TABLE IV: Number of events per bin versus analysis binning. 250,000 total detections are expected with the 
default experiment. 



ratio of at least 1:1 is desirable. The experiment may benefit from a second detector at longer baselines of 
10-20 m to improve sensitivity to the regions of small Am^. 




FIG. 17: Reactor and detector parameters relevant for covering the suggested parameter space. These graphs 
indicate the direction in which the sensitivity curve moved when reactor (left) and detector (right) parameters are 
improved or adjusted. 

Figure [17] qualitatively summarizes the impact of various reactor and detector parameters on the parame- 
ter space covered by very short baseline reactor X/g experiments. The relative important of these parameters 
is not weighted in this figure, and some experimental variables such as the precise core shape have only 
secondary effects on the experiment. 

• Total Statistics: The per-bin statistical uncertainty can be a limiting quantity when searching for 
oscillations over numerous energy and position bins. By utilizing a core with the highest possible power 
and a detector with the highest possible efficiency, cross-section, or proton density, one can maximize 
this quantity and considerably improve an experiment's sensitivity to short-baseline oscillations at all 
Am^ values. 

• Detector length: A large detector length increases an experiment's ability to resolve oscillations 
with position in addition to spectral distortions in energy. This effect can increase overall sensitivity 
at most Am^ values while extending the range of high sensitivity to lower values of Am^. 
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• Detector-reactor distance: The closest reactor-detector distance rmin determines the Am^ range 
of highest sensitivity. To achieve optimized sensitivity Vmin should be paired with a detector length d 
of equal or larger magnitude to allow sampling of a broad range of the oscillation period. In addition, 
statistics naturally increase as rmin is decreased. 

• Detector resolution: Oscillations at higher Am^ are only visible if resolutions and bin sizing are 
smaller than the oscillation itself. Sensitivity to values of Am^ approaching this threshold are signifi- 
cantly reduced. 

• Background: The S:B ratio is crucial for the success of the experiment. Small S:B ratios make 
it difficult to resolve oscillation effects above statistical background fluctuations and uncorrelated 
background uncertainties. For a given S:B, various background spectral shapes have similar impact 
on the experiment. However, precise knowledge of the backgrounds and their distribution in energy 
and position are critical for the experiment's sensitivity and for demonstrating the observation of 
neutrino oscillation. 

Other variables, such as core size, shape, and fuel type will not likely drive experimental design at most 
mass splittings. 



B. Energy Versus Energy+Position Measurements 



Two classes of detectors have been proposed for very short baseline measurements at reactors. The first 
consists of m'^-sized experiments that will look for oscillation-related distortions in the reactor energy 
spectrum. The second class is larger multi-m^ experiments with good energy and position resolution that are 
sensitive to both oscillations in position and energy. The sensitivity of these two types of experiments can 
be compared using the default detector-reactor arrangement for each case while varying the total detector 
length: the energy-only experiment is defined to have aim detector length, while the energy-|-position 
experiment is defined to have a 3 m detector length. A comparison of the sensitivities of three configurations 
is given in Figure |18| a 3 m detector with good position and energy resolution, aim detector with good 
position and energy resolution, and aim detector with no position reconstruction. 

Experiments with energy-|-position information provide significantly better sensitivity to oscillations than 
energy-only experiments. Observing distortions of the energy spectrum at multiple baselines will lessen the 
effect of uncertainties in the energy spectrum shape. If viewed at a single baseline, an energy spectrum 
distortion can be more easily described without oscillations by nuisance parameters for individual energy 
bins. This problem is amplified if oscillations occur at high energies, where the reactor spectral uncertainties 
are larger. In contrast, distortions from oscillation occur at different energy values for different baselines, an 
effect not easily neutralized with the comparatively rigid energy spectrum uncertainty nuisance parameters. 
Large differences also exist between experiments utilizing energy -|-position information in 1-m and 3-m long 
detectors. In addition to the difference in statistics, the 1 m detector samples a shorter portion of the 
oscillation period as discussed in Section |VB| Even better sensitivity is achieved by an experiment with 



two default detectors at distances of 4 m and 12 m as shown Figure 18 



IX. CONCLUSION 



New experiments at very short distances from reactors have the potential to make a precision measurement 
of the reactor antineutrino spectrum, resolve the "reactor anomaly" , and probe a large fraction of the 
sterile neutrino oscillation parameter space. This paper explores the experimental variables of a short- 
baseline reactor experiment and resulting experimental sensitivity. The reactor power, detector length, 
reactor-detector distance, energy and position resolutions of the detectors, and background to signal ratio 
are amongst the key parameters for such an experiment. 

Radially-extended detectors with good position and energy resolution provide a systematically robust 
approach to the search for neutrino oscillations at very short baselines. A highly-extended (>5m) detector or 
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FIG. 18: Comparison of sensitivities for oscillation experiments with energy-only or energy+position information. 
All experimental parameters, including closest reactor-detector distance, are identical, except the total detector 
length. Also pictured is an experiment utilizing an additional similar detector at a longer baseline of 12 m. This 
extra detector greatly increases the available range of baselines, significantly increasing the experimental sensitivity 
at lower Am^ values. 



multiple detectors at different baselines may expand the sensitivity of the experiment to all relevant regions 
of Am^. Furthermore, the use of moveable detectors would allow deployment at multiple radial locations, 
which may mitigate the possible effects of spatially varying backgrounds and will clearly demonstrate the 
effect of neutrino oscillations. 



Acknowledgments 



We thank Henry Band, Jeff Cherwinka, and Randy Johnson for useful discussions and comments. This 
work was prepared with support from the Department of Energy, Office of High Energy Physics, under grants 
DE-FG02-95ER40896 and DE-FG02-84ER-40153, the University of Wisconsin, and the Alfred P. Sloan 
Foundation. 



[1] T. Mueller et al, Phys. Rev. C 83, 054615 (2011), hep-ex/1101.2663. 

[2] A. Hayes, Reactor Antineutnno Flux and the Anomaly, AAP 2012, Honolulu, Hawaii, USA (2012), http: 

//www . phys . hawaii . edu/~ j gl/AAP/AAP_2012_Agenda . htm 
[3] P. Huber, Phys. Rev. C 84, 024617 (2011), hep-ex/1106.0687. 
[4] M. Fallot et al, Phys. Rev. Lett. 20, 202504 (2012), nucl-ex/1208.3877. 
[5] G. Mention et al, Phys. Rev. D 83, 073006 (2011), hep-ex/1101.2755. 
[6] T. Mueller et al, Phys. Rev. D 64, 112007 (2001). 
[7] A. Aguilar et al, hep-ex/1207.4809 (2012). 
[8] C. Giunti and M. Laveder, Phys. Rev. C 83, 065504 (2011). 
[9] E. Komatsu et al, Astrophys. J. Suppl. 192, 18 (2011), astro-ph/1001.4538. 
[10] F. An (Daya Bay), Phys. Rev. Lett. 108, 171803 (2012), hep-ex/1207.6632. 



22 



[11] Y. Abe (Double Chooz), Phys. Rev. Lett. 108, 131801 (2012), hep-ex/1207.6632. 
[12] J. Ahn (RENO), Phys. Rev. Lett. 108, 191802 (2012). 

[13] K.M. Heeger, Future of Reactor Experiments, Neutrino 2012, Kyoto, Japan (2012), | |http: //neu2012 .kek. jp/ 
neu2012/progrcLimne .html 



[14] T. Lasserre, Neutrino: Mass, Nature and Number, CIPANP 2012, St. Petersburg, Florida, USA (2012), [httpT 

//cipaiip2012 . triumf . ca/timet able .html 
[15] N. Bowden, New Short Baseline Reactor Experiments, SNAC 2011, Blacksburg, Virginia, USA (2012), [htt^ 

/ / www . cpe . vt . edu/ snac/pr ogr am . html 
[16] A. Serebrov et al, liep-pli/1205.2955 (2012). 

[17] V. Kopeikin et al. Atom. Phys. Nucl. 67, 1892 (2004), hep-ph/0410100. 
[18] Y. Abe (Double Chooz), hep-ex/1210.3748 (2012). 

[19] Oak Ridge National Laboratory, The High Flux Isotope Reactor at Oak Ridge National Laboratory (2012), 

http : //neutrons . ornl .gov/f acilities/HFIR 
[20] National Institute of Standards and Technology, NIST Center for Neutron Research (2012), |http: //www .ncnrT] 

nist .gov 

[21] Idaho National Laboratory, Advanced Test Reactor: Meeting US Nuclear Energy Research Challenges (2012), 

http : //www. ncnr .nist .gov 
[22] H. Kwon et al., Phys. Rev. D24, 1097 (1981). 

[23] San Onofre Nuclear Generating Station (2012), http : //www. songscommunity. com/default . asp?from=SONGS 
[24] V. Kopeikin et al, Phys. Atom. Nucl. 69, 185 (2006), hep-ph/0412044. 
[25] N. Bowden et al, ins-det/0612152 (2006). 
[26] T. Lasserre, Personal communication (2012). 

[27] Bryce Littlejohn, Ph.D. Thesis, University of Wisconsin - Madison (2012). 
[28] B. Ackhar et al, Nucl. Phys. B434, 503 (1995). 
[29] A. Afonin et al, JETP Lett. 54, 253 (1991). 
[30] A. Afonin et al, JETP Lett. 45, 247 (1987). 



[31] NIST Center for Neutron Research, 2012 Operations Schedule (2012), http://www.ncnr.nist.gov/ns. 



[32] Oak Ridge National Laboratory, HEIR Operating Forecast & Planning Schedule (2012), http : //neutrons . 

ornl . gov/ f acilit ies/HFIR/pdf /rollinglSmonth . pdf 
[33] ATR National Scientific Users Facility, FY 2009 Advanced Test Reactor National Scientific User Facility Users' 

Guide (2009), https : / / secure . inl . gov/atrproposal/documents/ATRUsersGuide .pdf 
[34] F. An (Daya Bay), hep-ex/0701029 (2007). ' ~ ' 

[35] F. An (Daya Bay), Nucl. Inst. Meth A685, 78 (2012), det-ins/1202.6181. 
[36] J. Beringer et al (Particle Data Group), Phys. Rev. D 86, 010001 (2012). 
[37] C. Giunti and M. Laveder, hep-ph/1109.4033 (2011). 



[38] San Onofre Nuclear Generating Station, Anatomy of a Nuclear Power Plant Planned Outage (2012), http: 

//www. songscommunity . com/docs/SONGS-Planned-Outage .pdf 
[39] Ross T. Ridcnoure, Docket No. 50-361 60-day Post Cycle-15 Refueling Outage Reactor Pressure Vessel Head 



Inspection Report for San Onofre Nuclear Generating Station, Unit 2 (2008), http://pbadupws.nrc.gov/docs/ 
ML0808/ML080850059 . pdf 

[40] D. Ridikas, G. Fioni, P. Goberis, O. Deruelle, M. Fadil, F. Marie, S. Rottger, On the Fuel Cycle and Neu- 



tron Fluxes of the High Flux Reactor at ILL Grenoble (2002), http://cdsweb.cern.ch/record/540488/files/ 
sis-2002-086.pdf 

[41] Alexander Glaser, Ph.D. Thesis, Technical University of Darmstadt (2005). 
[42] P. Vogel and J. Engel, Phys. Rev. D39, 3378 (1989). 
[43] K. Schreckenbach et al, Phys. Lett. B 160, 325 (1985). 
[44] A. Hahn et al, Phys. Lett. B 218, 365 (1989). 
[45] C. Lewis, personal communicaion (October 2012). 
[46] F. Boehm et al (Palo Verde), Phys. Rev. D64, 112001 (2001). 



